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Abstract 

Dynamical mass generation of a two-component fermion in QED% with a Chern- 
Simons term is investigated by solving the Schwinger-Dyson equation formulated 
in the lowest ladder approximation. Dependence of the dynamical fermion mass 
on a gauge- fixing parameter, a gauge coupling constant, and a topological mass is 
examined by approximated analytical and also numerical methods. The inclusion of 
the Chern-Simons term makes impossible to choose a peculiar gauge in which a wave 
function renormalization is absent. The numerical evaluation shows that the wave 
function renormalization is fairly close to 1 in the Landau gauge. It means that 
this gauge is still a specific gauge where the Ward-Takahashi identity is satisfied 
approximately. We also find that the dynamical mass is almost constant if the 
topological mass is larger than the coupling constant, while it decreases when the 
topological mass is comparable to or smaller than the coupling constant and tends 
to the value in QED3 without the Chern-Simons term. 
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I. Introduction 



In (2+l)-dimensions, generally in odd- dimensional space-time, a peculiar type of quantum 
field theory exists. That is a theory with a Chern-Simons term in its action. |]T], The 
Chern-Simons term may be added to the action by hand or induced in an effective action 
by a vacuum polarization effect of a fermion as the anomaly ||. Our understanding on 
the theory has been progressed especially through studies of the topological structure due 
to the Chern-Simons term. || 

On the other hand, two important phenomena were discovered in condensed matter 
systems. One is the quantum Hall effect || and the other is the high-T c superconduc- 
tivity 0. Both phenomena may be thought of as realization of macroscopic quantum 
effects in planar systems in which electrons have a strong correlation. It is challenging to 
ask whether and how the Chern-Simons term is related to these phenomena discovered in 
the two-dimensional electron systems. These situations have activated the study of the 
(2+l)-dimensional quantum field theories. 

In (2+l)-dimensional gauge theories, the behaviour in the infrared region seems to be 
unstable at least in the perturbative treatment. The Chern-Simons term gives the gauge 
field a topological mass without breaking gauge symmetry so that the term may rescue 
the theory from the infrared catastrophe. This is original motivation of including the 
Chern-Simons term in the action. Reminding those studies, we think that one of the 
most important problems which is still not clear is what is a role of the Chern-Simons 
term in a nonperturbative region. It is our purpose to clarify how the Chern-Simons term 
affects nonperturbative dynamics. 

The dynamical mass generation of the fermion is an important phenomenon induced 
by the nonperturbative effects. In even- dimensional space-time, chiral symmetry forbids 
the mass of the fermion. The nonperturbative effects can break this symmetry giving the 
mass to the fermion. As is well known, we cannot define the chiral symmetry for the 
two-component fermion which belongs to the irreducible spinor representation in (2+1)- 
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dimensions. Instead parity symmetry (P 2 ) forbids the two-component fermion having 
the mass. On the other hand, for the /owr-component fermion which is composed of 
two two-component fermions, we can define a kind of "chiral" symmetry. The "crural" 
transformation is defined by a combination of the parity and Z 2 flavour transformations. 
(P A = P® Z 2 ) 

Many studies on the dynamical mass generation in QED3 have been done already 
by using the Schwinger- Dyson technique [H] or the lattice Monte-Carlo method ||. In 
almost of previous works, these authors have studied the dynamical mass generation of 
the four-component fermions in the theory without the Chern-Simons term, which is a 
spontaneous symmetry breaking of P4 but not P 2 . These studies have been extended 
to the case with the Chern-Simons term but the fermion considered there has been the 
four-component one. [|TIJ 

The theory with iV four-component fermions is equivalent to the theory with 2N 
two-component fermions. In this sense, the theory with the four-component fermions 
can describe only the case with even number of flavours (counting the two-component 
fermions). It is interesting to ask what happens in the theory with odd number of flavours. 
This also is another motivation of the present work. The simplest version of the theory 
with odd number of flavours is the case with a single flavour. The first work concerned with 
the dynamical mass generation of the single two-component fermion in QED3 without the 
Chern-Simons term was done in Ref . |IT] . The symmetry which was studied in the work 
is P 2 but not P4. In this paper, we extend the work to the case with the Chern-Simons 
term. 

Thus we consider the dynamical mass generation of the single two-component fermion 
in QED 3 with the Chern-Simons term. The Schwinger-Dyson equations are formulated 
in the lowest ladder approximation, which are two coupled integral equations. Due to 
the Chern-Simons term, a drastic change happens in the structure of the Schwinger- 
Dyson equations. While the wave function renormalization is absent in the Landau gauge 
in QED% without the Chern-Simons term as far as we consider in the lowest ladder 
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approximation, the inclusion of the Chern-Simons term makes impossible to choose a 
peculiar gauge in which the wave function is not renormalized. Because there are no ad 
hoc reasons that the Landau gauge is specific in QED3 with the Chern-Simons term, we 
should study the coupled integral equations for various values of the gauge parameter. 

Our strategy in this work are as follows: First we solve the coupled integral equations in 
an approximated analytical method which is very crude but may be available for reference 
in a more complete numerical calculation. After that, the coupled integral equations are 
solved numerically. In these analyses, we study dependence of the dynamical mass on the 
gauge-fixing parameter, the coupling constant, and the topological mass. 

This paper is organized as follows: In Sec. II, we explain the model which we use. 
The Schwinger-Dyson equations in the lowest ladder approximation are derived in Sec. 
III. In Sec. IV, the equations are analyzed by the approximated analytical method. The 
numerical analysis of the equations are presented in Sec. V. We conclude our results with 
discussions in Sec. VI. In Appendices A and B, some useful formulae are summarized. 



II. Maxwell-Chern-Simons QED3 

The model which we consider is an extended version of the usual three-dimensional quan- 
tum electrodynamics (QED3). It has both of the Maxwell and Chern-Simons terms as the 
action of the U(l) gauge field . We call this theory the Maxwell-Chern-Simons QED 3 . 
The gauge field is coupled to the two-component Dirac fermion. The Lagrangian density 
of the theory is given by 

£ = --F^ + ^"A^A,, - ^{d^f +$(i@-e f)^ . (1) 
4 2 2a 

The second term in the right-hand side of Eq.(|l]) is the so-called Chern-Simons term. It 
is well-known that the term gives the gauge field the mass fi without breaking the gauge 
symmetry, fi is called the topological mass because the Chern-Simons term has topological 
meaning as a secondary characteristic class M, P2fl . a is the gauge- fixing parameter and 



e is the gauge coupling constant, ip is the two-component fermion field which belongs to 



the irreducible spinor representation in (2+l)-dimensions. The Dirac matrices are defined 
by 7° = 03, 7 1 = icrxj'j 2 = ia 2 with diag(g^ u ) = (1, —1, —1) where <Tj's (i=l, 2, 3) are the 
Pauli matrices. The 7 M 's satisfy relations as {7^, Y} — 2g ,M ", 7^7" = —ie' MUP/ y p + g pv and 
tr[7 M 7^] = 2g fJ,u . In this representation, there does not exist a matrix which anti-co minutes 
with all of 7 M 's so that we cannot define the chiral transformation. This is a specific aspect 
of the odd-dimensional space-time. In even-dimensions, the chiral symmetry requires 
that a fermion is massless. In odd- dimensions, the chiral symmetry itself does not exist. 
Instead, the mass term of the fermion is forbidden by parity symmetry. The parity 
transformation is defined as 

x = (t,x,y) -> x' = (t,-x,y) , 

ip(x) -> 7V<y) , 

A°(x) -+ A°(x') , (2) 
A 1 (x) -> -A\x') , 
A 2 (x) -> A 2 (x') . 

Under the parity transformation, the mass term of the fermion and the Chern-Simons 
term change their signs. Thus the mass terms of both the fermion and the gauge field are 
forbidden by the parity symmetry. We study how the breaking of parity by the topological 
mass affects the mass generation of the fermion. 

III. Schwinger-Dyson Equations 

As a non-perturbative method to evaluate the dynamical mass of the fermion, we use 
here the Schwinger-Dyson technique. The Schwinger-Dyson equation for the fermion 
self-energy E(p) is written as 

E(p) = (^e) 2 J^'fiS , F (k)r'(k,p-k)iD' ltv ( P -k) . (3) 
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T u (k, p — k) is a full vertex function and D'^ u (p — k) is a full propagator of the gauge field. 
S' F is the full propagator of the fermion field which is written as 

where A(p) and B(p) are functions of y/p^p 11 in taking into account the relativistic nature, 
while S(p) depends on p^s. A(p)~ l is the wave function renormalization and B(p)/A(p) 
is a mass induced by dynamical effects at the momentum scale p. The so-called dynamical 
mass m phys is defined by m phys = B(0)/A(Q) as usual. It is useful to notice the relations 



as 



tr[E(p)) = -2iB(p) , tr[jfE(p)] = 2i(A(p) - l)p 2 . (5) 

To proceed our analysis of Eq. (Q) further, we need to introduce any suitable approxi- 
mation. In this paper, we limit ourselves to use the lowest ladder approximation where the 
full propagator of the gauge field and the full vertex are replaced by the free propagator 
and the tree vertex respectively as 

iD'^{p-k)^iD^{p-k) , T\k,p-k)^Y ■ (6) 

The analysis beyond the ladder approximation will appear elsewhere. The free propagator 
of the gauge field is derived from the Lagrangian density (|l|) as 

iDr(p) = [gT - ^f] + »-^- 2 -^ VP V P - ia^f ■ (7) 

p 2 — fi 2 \ p 2 J p 2 — \i 2 p 2 p A 

Thus the Schwinger-Dyson equation in the lowest ladder approximation becomes 

S(p) = He) 2 J ^l»iS' F {k)l V iD» v (p-k) . (8) 

We substitute Eqs.(f|) and (0) into Eq.® and use Eq.® in producing two coupled 
equations. After taking the traces, we obtain the coupled integral equations as 

p 2 — k 2 — (p — k) 2 



, * ie 2 r d 3 k I 

B \p) = — 7T 



2 J (2ir) 3 A(k) 2 k 2 - B(k) 2 



-2fiA(k) 



{(p - k) 2 - // 2 }(p - kf 
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+ AB(k) 



A{p) 



1 - 



le 



(p — k) 2 — p' 
2 - d 3 k 



+ 2aB{k) 



(p - k)< 



(9) 



2p 2 J (2tt) 3 A(k) 2 k 2 - B(k) 2 
2pB(k) 



A(k) 



(P 2 



„2\2 



(p — k) 2 — p 2 [ (p — ky 



-(p-kf 



k 2 



(p — k) 2 — p 2 \ {p — ky 



+ 1 



aA(k) 



(p 2 — k 



2\2 



p 2 + k 2 



(p - ky (p - ky 



(10) 



Now we change the metric to the Euclidean one by the Wick rotation as (k°, k) —>■ 
(ik°, k) and {p°,p) — > (ip°,p). Then k 2 and p 2 are replaced by — k 2 = —(k ) 2 — (k 1 ) 2 and 
— p 2 = —(p ) 2 — (p 1 ) 2 - After that, we transform the integral variables /c M 's to the polar 
coordinates (k,9,(f>). The angular integration on 9 and <fi can be done explicitly. (The 
useful integral formulae are listed in Appendix A.) Finally we obtain the coupled integral 
equations which contain only the integration on the radial variable k, 



B{p) 



^-K 2 p JO 



dk 



k 



A(k) 2 k 2 + B{k) 2 



aB{k) - -{p 2 - k 2 )A(k)\\n 



(p + ky 
(p — ky 



+ 



A{p) 



1 + 



-Go 2 - k 2 )A{k) + fiA(k) + 2B(k) In 



(p + k) 2 + n 2 
(p — k) 2 + jj? 



Ill 



]ir 2 p 3 Jo 



dk 



k 



A(k) 2 k 2 + B(k) 2 



1 



-2pk(a + l)A{k) + { ^- 2 {p 2 - k 2 ) 2 A{k) 



+ 



1 (p 2 - k 2 )B(k) + \a(j? + k 2 )A{k)} In |±|^ + {^ 2 A(fc) 



(p 2 - k 2 ) 2 A(k) + ^(fc) - i(p 2 - A: 2 )!^)} In |^|^ 



;i2) 



In the successive Sees. IV and V, we solve these equations by an approximated analytical 
method and also numerically by using an iteration method. 

Notice that the parameters e 2 and p have the dimensions of mass. We may rewrite 
Eqs. (|TT|) and flT2" ) to dimensionless forms by defining jx = /i/e 2 , b(x) = B(e 2 x)/e 2 and 
a(x) = v4(e 2 a;) where x is a dimensionless variable defined by p = e 2 x. The equations are 
just the ones obtained by setting e 2 = 1 and replacing (A, B, p) to (a, b, p) in Eqs. (|TT|) 
and (|12"D. The theory is controlled by only one dimensionless parameter p. After solving 
the dimensionless equations for a(x) and b(x), we can convert them to A(p) and B(p). 
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IV. Approximated Analytical Studies 
IV- A. fi limit 



We can check easily that Eqs. (|TT|) and (0) reduce to the Schwinger- Dyson equations in 
QED3 without the Chern-Simons term if we put the topological mass fi equal to zero. In 
fact, taking the limit as /i — > in Eqs. ( |TTD and (0), we obtain 



B(p) 
A(p) 



(a + 2] 



8ir 2 p Jo 

,2 



dk 



kB{k) 



■In 



(p + A;) 5 



A(/c) 2 /c 2 + B(k) 2 (p — A;) 2 ' 
jb4(Jfe) 



a 47rVio A(k) 2 k 2 + B{k) 2 



p 2 + k 2 1 (p + A;) 2 
^"^4 n (p~^ 



(13) 
(14) 



which are the Schwinger- Dyson equations in the lowest ladder approximation derived in 
QED 3 without the Chern-Simons term. 2 We can see that there exists the specific gauge 
where the wave function renormalization is absent. Thus in the Landau gauge (a = 0), 
Eq. ( PI) gives us the simple solution as A{p) = 1 and the problem reduces to solve Eq. 
([13]) with A(p) = l. 

In the case with the Chern-Simons term, as is seen in Eqs. (|Tl~l) and (|i~2|), there does 
not exist such a specific gauge where the wave function is not renormalized. As far as 
we cannot find a self-evident reason that the Landau is still specific in QED 3 with the 
Chern-Simons term, it must be fair to study Eqs. (]TJ) and flT^D for various values of the 
gauge-fixing parameter a. 



IV-B. Small p Expansion and Perturbative Results 



At first sight, Eqs. (11) and fll2[) seem to be dangerous in the limit p — * 0. But we can 
check that this limit is well-defined. In the region p « 0, Eqs. ([11]) and ( 12 ) are written as 



B(p) 



A(p) 



IT 2 JO 



dk 



k 



A(k) 2 k 2 + B(ky 



kfI A(k)+(^ + 



k 



k 2 + fi 



+ 0(p 2 



—T / dk 



ir 2 Jo A(k) 2 k 2 + B(ky 



1 J a 

3\k~ {k 2 + l i 2 ) 2 



\2k k 2 + /j, 2 t 
2k^ 2 ^ 



B{k) 



(15) 



A{k) 



In there were some typing mistakes. We give the correct equations in Eqs. (|13[) and (|14|). 



Sk (fc z + /i z j z 



(See Appendix B.) From Eqs. (|T3D and ([16|), we can derive the result which is obtained 
in the lowest-order of perturbation. By setting A(k) = 1 and B(k) = and performing 
the integration on k, we have 

B(0) = f M , 40) = 1 - £M + ^ , (17) 

27T /i D7T /x 2 37^ e 

where e is the infrared cutoff in the integration on k. It should be noticed that B(0) 
depends on only the sign of fi. This also may be a specific aspect in the (2+l)-dimensions. 
The dependence of A(0) on \x shows that only the Landau gauge is free from the infrared 
divergence. On the other hand, A(0) is singular at \x = so that the theory with the 
Chern-Simons term may not be smoothly connected to the theory without the Chern- 
Simons term in the lowest-order perturbation. 

IV-C. Constant Approximation 

Before proceeding to a numerical analysis, it is very useful if we can estimate A(0) and 
B(0) analytically even under a fairly crude approximation. The kernels of these integral 
equations are dumped rapidly as the integral variable k increases so that the contribution 
from k ~ is the most dominant one in the integrals. We approximate A(k) and B(k) by 
A(0) and B(0) in the integrals. We call this approximation "the constant approximation" . 
Of course this approximation might be too crude for our purpose and we only use the 
result as reference in the numerical analysis. Under this approximation, we can perform 
the remaining radial integration and have the simple algebraic equations as 

e 2 1 a 

= T,W) [ 2 +l) • (18) 

- ' (19) 

where we have considered the case of /i > 0. (See Appendix B for the details.) Solving 
the coupled algebraic equations ( ^8|) and (|T9|) , we obtain 



e 2 e 2 



B ^ = 2-n + U; a ' (20) 
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A(0) = 1 + -^ • (21) 
a + o 

From Eqs. (|20|) and (|21|), we can see that the dependence of -B(O) and A(0) on the 
gauge-fixing parameter, the coupling constant and the topological mass has the following 
peculiar features: 

1) Dependence on the gauge- fixing parameter 

B(0) depends linearly on a. It is suggestive that A(0) is singular at a = — 6 where 
B(0) vanishes. In the Landau gauge (a = 0), A(0) = 1 and B(0) = e 2 /2ir. A(0) = 1 
is favourable for us because A(p) = 1 means that the Ward-Takahashi identity is 
satisfied. 

2) Dependence on the coupling constant 

A(0) does not depend on e 2 . It means that the deviation of A(0) from 1 is inde- 
pendent of the coupling constant. This is crucially different from the perturbative 
result given by Eq. (|17D where the deviation is proportional to e 2 . On the other 
hand, B(0) is proportional to e 2 . 

3) Dependence on the topological mass 

We recognize that Eqs. (|20|) and ( |2~ll) are independent of the topological mass \i. In 
fact, if we apply the constant approximation to the case without the Chern-Simons 
term, we obtain the same results as Eqs. (|20D and (pip . It means that the amount 
of the explicit parity breaking in the gauge sector by the topological mass does not 
affect the dynamical mass in the fermion sector. 

Now we proceed to a more precise numerical evaluation in the next section. 

V. Numerical Analysis 
V-A. Non-trivial Solutions 

We solve the two coupled integral equations ( PD and (^) numerically by using a method 
of iteration. First we substitute trial functions into A(k) and B(k) in the right-hand 
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sides of Eqs. (|TT|) and (|T2"D and then calculate the integrals numerically. The outputs so 
obtained, A(p) and B(p), are substituted back to the right-hand sides until the outputs 
coincide with the inputs. Finally we obtain convergent functions A(p) and B(p), which 
satisfies the integral equations, if there exist any solutions in Eqs. (|ll|) and ([12|). 

We obtain the non-trivial solutions for the various values of the gauge parameter a. 
For e = 1.0 and fi = 1.0, A(p) J s in the cases of a = 0, 1,2,3 are shown in Fig. 1. The 
corresponding B(p)'s are presented in Fig. 2. We find a very interesting feature that A(p) 
is fairly close to 1 in the Landau gauge (a = 0). 

In the case of QED 3 without the Chern-Simons term, A(p) is exactly equal to 1 in 
the Landau gauge under the lowest ladder approximation [|11|] However, in the case of 
QED3 with the Chern-Simons, there may be no apparent reason that A(p) = 1 in the 
Landau gauge. It is surprising that the numerical calculation of so complicated integral 
equations results A(p) « 1 in the Landau gauge. There might be a simple reason of 
explaining a peculiarity of the Landau gauge. 

V-B. Dependence on the gauge parameter 

The dependence of A(0) and B(0) on the gauge-fixing parameter a is shown in Figs. 3 
and 4. In the region a > 0, the numerical result is consistent with the one obtained in 
the constant approximation. 

For a < 0, the numerical iteration procedure does not converge, but cycles between 
two or more functions, none of which are solutions to Eqs. (|11|) and (|12"D . This appears to 
be a manifestation of the well known "doubling route to chaos" frequently exhibited by 
non-linear iterative algorithms. Unfortunately this prevents us from obtaining numerical 
solutions for negative values of the gauge parameter. 

A plot of the gauge invariant condensate < ipip > as a function of a is helpful as a way 
of indicating to what extent gauge symmetry is broken by the bare vertex approximation. 
The condensate is defined by < ipip >= —i\im x ^ trS' F (x) where S' F (x) is a propagator 
3 Outside the lowest ladder approximation, it is known that A(p) differs from one in the infrared, pi 



10 



in real space-time coordinates. Using the Fourier transformation and the Wick rotation, 
we obtain 



On the other hand, a position of the pole of the Minkowski propagator also is gauge 
invariant. To know the position, we have to perform an analytic continuation from our 
Euclidean results to the Minkowski ones. But it is a hard task because what we know in the 
Euclidean analysis is numerical. In fact, it is difficult to derive a full analytic properties 
from the numerical data. Instead, B(0)/A(0) is usually considered as an approximated 
value of p 2 at the pole. We have shown the a-dependence of the dynamical fermion mass 
mphys — B(0)/A(Q) and also < ipifj > in Fig. 5. It shows that the a-dependence may be 
considered to be fairly weak, compared with the result in the constant approximation. 

The results obtained in Secs.V-A and V-B suggest that the Landau gauge is still the 
best gauge. Hereafter we present mainly the results obtained in the Landau gauge. 

V-C. Dependence on the topological mass and the coupling con- 
stant 

Now we investigate how the dynamical fermion mass depends on the topological mass and 
the coupling constant. What we are most interested in is the dependence of the dynamical 
fermion mass on the topological mass of the gauge field. In the constant approximation, 
it has been shown that both A(0) and B(0) do not depend on the topological mass. Is 
this true in the more precise numerical evaluation? 

In Fig. 6 we show the dependence of a(0) on the dimensionless parameter fi which 
is defined in Sec. III. We can see that the deviation of a(0) from 1 is less than 1 %. 
We may say that a(0) is almost equal to 1 in all region of jl. It means that the fi- and 
e 2 -dependence of A(0) is extremely weak. 

On the other hand, the /i-dependence of 6(0) is presented in Fig. 7. It show that 
though B(0) is almost constant in the region of ji > e 2 , it decreases as fi does if \i is 




(22) 
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comparable to or smaller than e 2 . The e 2 -dependence of -B(O) is less explicit in Fig. 7 so 
that we show the explicit e 2 -dependence of B(0) in Fig. 8. B(0) increases as the coupling 
constant becomes larger. The e 2 -dependence of B(0) is almost linear when e 2 is smaller 
than /i. This is consistent with the constant approximation. But the slope is smaller than 
the one obtained in the constant approximation and the line curves downward slightly as 
the coupling constant becomes larger, which may be a non-perturbative effect. 

It should be noticed that for very small or large values of the topological mass com- 
pared with the coupling constant, there appears a technical difficulty in the numerical 
calculation as follows: In principle, we do not need to cut off the region of the energy- 
momentum integrals in Eqs. (|TT|) and (|T2|). But in the practical prescription of the 
numerical integration, the infinite range of the integration must be replaced by a finite 
one so that it is not avoidable to introduce cut-off parameters. The cut-off dependence of 
the results must be checked carefully. When the topological mass fi takes very small or 
large values compared with e 2 , it has been found that the integration region must be taken 
wider enough to get reliable results. It needs much machine power. In getting our results, 
we have checked the absence of the cut-off dependence in detail in the parameter region 
adopted here (/i/e 2 = 10~ 2 ~ 10 4 ). Therefore we conclude that the behaviour of A(0) 
and B(0) mentioned above, especially the decreasing of B(Q) in the region of /i/e 2 < 1, is 
not due to the cut-off. 

VI. Conclusion and Discussion 

We have investigated the dynamical mass generation of the single two-component fermion 
in QED 3 with the Chern-Simons term. The coupled Schwinger-Dyson equations for A(p) 
and B(p), where the full fermion propagator is defined by S' F (p)~ 1 = A{p) p 1 — B(p), have 
been formulated in the lowest ladder approximation. When the Chern-Simons term is 
included, there does not exist the specific gauge where A(p) is automatically equal to 1 as 
the case without the Chern-Simons term. We examine the dependence of the dynamical 
mass on the gauge-fixing parameter a, the coupling constant e and the topological mass 
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fj, by using the approximated analytical and also the numerical methods. 

The coupled integral equations have been solved first analytically in the constant 
approximation where A(p) and B(p) have been replaced by A(0) and B(0) respectively. 
We have found that -B(O) is proportional to a and also to e 2 . On the other hand, A(0) is 
independent of e and singular at a = —6 where -B(O) vanishes. We also have found that 
A(0) and B(0) do not depend on the topological mass. 

Keeping these facts in mind we have proceeded to solve the set of the integral equations 
by using the numerical method. The dependence of A(0) and -B(O) on a and e are almost 
consistent with the results derived in the constant approximation. -B(O) depends on a 
linearly. The e 2 -dependence of -B(O) also is almost linear but the value of the slope becomes 
smaller than the one in the constant approximation as e 2 increases. The e-dependence of 
A(0) is almost absent. 

To check to what extent gauge symmetry is spoilt by the bare vertex approximation, 
we have evaluated the gauge invariant condensate < xpip > and also B(0)/A(0). The 
result has shown that a-dependence of them may be considered to be fairly weak. 

Further we have discovered some novel features in the numerical analysis. First we 
have found A(p) « 1 in the Landau gauge. It is not entirely obvious why the numerical 
evaluation of so complicated integral equations results that ^4(0) is almost equal to 1 at a = 
0. We may expect that there is any simple reason. It also should be noticed that the trivial 



A(p) in the Landau gauge is obtained merely in the lowest ladder approximation. [I3| More 
sophisticated approximation will be considered in subsequent works. 

Secondly we have found the strange behaviour of A(0) and -B(O) in the region a < 
0. There has appeared the doubling signal in the iterations. The signal is a popular 
phenomena in a non-linear system. We wonder that the signal appeared here would have 
any physical or mathematical meaning or not. 

As the third, it seems that the dynamical mass of the fermion is almost constant if fi is 
larger than e 2 , while it decreases when /i is comparable to or smaller than e 2 . Why it is? 
In the case of fi ^ 0, we have found just one solution in our numerical analysis. On the 



13 



other hand, we know that there exist two solutions in the case of \x = [11|; one is trivial 
(B = 0) and the other one nontrivial (B ^ 0). We have extrapolated our numerical data 



B(0)/A(0) to fj, = by using the method of least squares [[U| and obtained Fig. 9. The 
numerical evaluations in QED% without the Chern-Simons term results B(0) = 0.104755 
(a = and A(0) = 1). The value obtained by the extrapolation is B(0)/A(0) = 0.105255 
so that both results may be consistent with each other. Thus, the solution in the case 
of fj, seems to approach the non-trivial one as \i —>■ 0. Of course, we should be 
careful before accepting this result. As mentioned in V-C, the numerical evaluation in 
the region of fi/e 2 << 1 is very difficult technically and the behaviour of the solution in 
taking the limit /x — > is still not clear. There might appear any critical behaviour such 
as a bifurcation of the solution. The best way to confirm the result is to give a proof 
by analytical method but it is very difficult because the integral equations are highly 
non-linear. We are now extending the parameter region to search any more non-trivial 
structure especially in the region of /i/e 2 << 1. The result will be reported separately. 
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Appendix A: Angular Integration 



We present here the formulae for the angular integration, which are used to derive Eqs. ([11]) 
and flT2|) in Sec. III. After the Wick rotation, Eqs. (|S]) and (P~D|) can be rewritten in the 
polar coordinates as 

^2 poo h% 



B(p) 



A{p) 



dk- 



8n 2 Jo A(k) 2 k 2 + B(k) 2 
+ 2^A{k)h + AB(k)h + 2aB(k)I 2 



2fi(p 2 - k 2 )A(k)I 



dk- 



k 2 



A(k){(p 2 -k 2 ) 2 I -h} 



8n 2 p 2 Jo A(k) 2 k 2 + B(k) 2 
+ 2fiB(k) {{p 2 - k 2 )I + h} - aA{k) {(p 2 - k 2 ) 2 h - (p 2 + k 2 )I 2 ) 



Each integrals J n 's (n=0, 1, 2, 3, 4) are calculated as 

1 



h = / 
Jo 



{(p - k) 2 + /i 2 }(p - k) 2 2/j 2 pk (p- k) 2 {(p + k) 2 + v 2 } 
1 



h 
h 
h 



sin 6d9 
sin 6d6 
sin Odd 



1 ln (p + k) 2 {(p-k) 2 + fx 2 } 



1 (p + A;) 2 + /i 2 



(p — A;) 2 + /i 2 2pA; (p — k) 2 + /i 2 
1 



1 ln ( P + k) 2 



o (p — k) 2 2pk (p — k) 2 

sin 6d6 



o 



{P ~ k)2 1 ^ P+ k) 2 -(p-k) 2 -^J P + k)2 + fi2 



(p — k) 2 + /i 2 2pA; 
1 If 1 



(p — k) 2 + /i 2 



1 



„ sm 9(19 (p _ fc )4 2 pA: | (p + A;) 2 (p - A;) 2 



Appendix B: p — ► limit and Constant Approximation 



The formulae which are used in Sec. IV are summarized here. The (p — > 0) limit is taken 
by the expansion formulae as 



In 



iP + k)2 ~ 4? + l(?)3 + 0(p 5 ) 



In 



(p — k) 2 
{p + k) 2 + fx 2 



k Tk' 

Akp 4 (A; 2 - 3 / u 2 )A; 



+ 



p 3 + 0(p 5 



(p-k) 2 + fi 2 k 2 + fi 2 3 (A; 2 + /i 2 ) 3 
In the constant approximation of Eqs. (|11|) and (|T2l), we replace the unknown functions 
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A(k) and B(k) to A(0) and B(0). Setting p = 0, we obtain 

e 2 r 
5(0) = € - 2 



IT' 



M(o) + (^ + i)5(o) [> Jo - /i 2 (M(o) + 5(0)) Jx 



A(0) = 1 + 4 



3 

|A(0) J - |(2M(0) + 5(0)) Ji + ^-(2M(0) + 25(0)) J 2 



where Jj's (i = 0, 1, 2) are given by 



J 



Jo 



o A(0) 2 k 2 + 5(0) 2 2,4(0)5(0) ' 

T f°° dk 1 7r 1 1 

Ji = / 



J 



2 



A(0) 2 fc 2 + £(0) 2 fc 2 + // 2 25(0) \^\A(0) + 5(0) ' 

rfA; 1 TT 1 2|/x|A(0) + 5(0) 



/o A(0) 2 A;2 + 5(0)2(A;2 + / i2)2 45(0) |/i| 3 (|/i|A(0) + 5(0)) 2 
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Figure 1: The nontrivial solutions of A(p) at e = 1.0 and /i = 
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Figure 2: The nontrivial solutions of B{p) at e = 1.0 and ji - 
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Figure 3: The a-dependence of A(0) at e = 1.0 and ji = 
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Figure 4: The a-dependence of B(0) at e = 1.0 and \x 
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Figure 5: The a-dependence of B(0)/A(0) and < tpip > at e = 1.0 and /i — 1.0 (0 < a < 
10) 



alpha=0 and e=1.0 



Numerical Result a(0) 
Constant Approximation 



--—.^^coe^e 



1 10 

mu/e A 2 



Figure 6: The /x-dependence of a(0) at a = and e = 1.0 
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Figure 7: The /i-dependence of 6(0) at a = and e = 1.0 
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Figure 8: The e 2 -dependence of B(0) at a = and fi = 1.0 (1.0 x 10 4 < e 2 < 
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Figure 9: The extrapolation of B(0)/A(0) in the region of fj, — 0.01 ~ 0.1 to /i — (a — 
and e = 1.0). 
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